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f~>) ' We study Friedmann-Robertson- Walker cosmological models with matter content composed of 

\ two perfect fluids pi and p2, with barotropic pressure densities pi /pi = cji = const and pijp^ = 

UJ2 = const, where one of the energy densities is given by pi — Cia°' + C2a^ , with Ci, C2, at and 
/3 taking constant values. We solve the field equations by using the conservation equation without 
breaking it into two interacting parts with the help of a coupling interacting term Q. Nevertheless, 
with the found solution may be associated an interacting term Q, and then a number of cosmological 
00 ' interacting models studied in the literature correspond to particular cases of our cosmological model. 

Specifically those models having constant coupling parameters a, /3 and interacting terms given by 
Q = oiHp^^j, Q = aHp^^, Q = aH{p^^, +Pde) and Q = 5iHp^^, +pHp^^, where p^M and p^E 
are the energy densities of dark matter and dark energy respectively. The studied set of solutions 
contains a class of cosmological models presenting a scaling behavior at early and at late times. On 
I the other hand the two-fiuid cosmological models considered in this paper also permit a three fiuid 

interpretation which is also discussed. In this reinterpretation, for fiat Friedmann-Robertson- Walker 
cosmologies, the requirement of positivity of energy densities of the dark matter and dark energy 
components allows the state parameter of dark energy to be in the range —1.37 < ^^de < —1/3. 
^ , Keywords: dark energy theory, dark matter theory 

I. INTRODUCTION 

> 

ly^ ' Recent observational data indicate that our Universe is undergoing accelerated expansion. The standard modern 
cosmological models consider the total energy density of the Universe to be dominated today by the densities of two 

■ components: dark matter (which has an attractive gravitational effect like usual matter), and dark energy (a kind of 
[ vacuum energy with a negative pressure), which drives the accelerated expansion [jy]. Even more, observational data 

fSj ■ seem to indicate that the Universe today may be dominated by an exotic kind of dark energy which has a very strong 
' negative pressure, denominated phantom fluid since it violates all energy conditions ^^l. The real nature of the dark 
. sector remains unknown. 

Most considered dark energy models present an accelerated expansion due to the presence of a quintessence (de- 
scribed by a canonical scalar field), or a phantom field (described by a scalar field with a negative kinetic term), 
among others. It appears in this kind of models that the energy density of dark energy is almost equivalent to that of 
^ . , the matter in the current or in recent times although they scale independently during all the cosmic evolution. This 

■ situation is referred to as the "coincidence problem" . To solve the coincidence problem many attempts have been 
done. Of special interest are interacting cosmological models since they may alleviate or even solve the coincidence 
problem fS]. We shall refer to this type of cosmological models below in this section. 

In the framework of General Relativity, for modeling the present state of the Universe, usually the cosmological 
models consider two cosmic fluids as sources for the Einstein field equations: one fluid for the dark matter sector, 
where is included the standard visible matter, and another cosmic fluid for the dark energy sector. Usually the dark 
matter is described as a pressureless ideal fluid, while the dark energy may be described by a cosmological constant, 
or perfect fluid or scalar fields 
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The consideration of two fluids in Einstein field equations leads to 

^A"' - ^Rgt^i^ = 1^ [Tlu + Tl^) : (1) 

where k, = SttG, and T^^ and T^^^ are the energy-momentum tensors of the two fluids. In a homogeneous and isotropic 
FRW universe 

ds^ = dt^ - aitf (^Y^^ - r''{de^ + sin^ddifi^)^ , (2) 
filled with two fluids and , the Friedmann equation is given by 

3jy2 = «(p,+pj- A. (3) 

From Eq. ([T]) it follows that both fluids together form a system that has a conserved four-momentum, implying that 
the two components and satisfy the following conservation equation: 

Pi+P2+3i?(Pi+P, +Pi+pJ = 0. (4) 

It is clear that the single two-fluid conservation equation (U) implies that the sum of the two fluids is conserved. In 
order to find solutions one may separate this single two-fluid conservation equation into two conservation equations 
for a system of two single fluids: 

p,+3H{p^+pJ=0. (5) 

The physical interpretation is that each of two fluid components is conserved, evolving separately according to standard 
conservation laws. However, this splitting into two conservation equations is a condition imposed on the conservation 
equation (|4]) and thus leads to a particular set of solutions to Einstein field equations. If we apply this picture to the 
dark sector we conclude that dark matter and dark energy are conserved separately, and if the pressures are given by 

p^{t) ^ uj^p^{t),p^{t) ^ uj^p^{t), (6) 

where and are constant parameters, then the solutions of Eqs. ([5|) are given by 

pi = pioa"'^'+"^',P2 = P2oa"'^'+"^', (7) 

implying that energy density of dark matter has the form p^^j^^ — Pdmo'^ 

However, since the real natures of dark matter and dark energy remain unknown, one can consider more general 
scenarios where dark matter and dark energy should not conserve separately and are coupled to each other. One 
coupling mechanism can be formally introduced into the Friedmann equations by defining an interacting term Q{t) 
in the following form: 

p, +3F(p, +pj =Q(i), 
p,+3H{p,+p,)^-Q{t). (8) 

Note that the case Q > is interpreted as a transfer of energy from fluid p^ to fluid p^ . Then, for the case Q < 0, we 
should have an energy transfer from fluid p^ to fluid p^ . In order to find solutions the form of this phenomenological 
interacting term Q must be postulated. Usually one of the following types is considered: Q = aHp^^^j , Q = aHp^^^ , 
Q = aHp^f^j -f fSHp^j^, Q ~ aH(p^j^j + p^^), where p^j,, and p^^ are the energy densities of dark matter and dark 
energy respectively. 

As we stated above, the introduction of this type of coupling between components of the dark sector offers to the 
standard cosmology a potential solution to the cosmic coincidence problem by requiring the ratio of matter and dark 
energy densities to be stable against perturbations at late times On the other hand, it also must be noticed that 
this type of coupling is motivated by considerations of high energy particle physics @ . 

It is clear that the assumption of a gravitational coupling between these two fluids in the form of Eqs. (HI is too close 
to the physical character of the two-fluid conservation equation since the total energy is conserved always satisfying 
the conservation equation Nevertheless, one can look for a general solution to Einstein equations by using the 
general properties encoded into the two-fluid conservation equation by choosing another more natural way to solve 
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the Friedmann equations ([3]) and (|4]), namely by postulating the form of one of the energy densities based on physical 
considerations and then find the whole solution by means of the conservation equation Q . In this paper this method 
will be applied. We shall consider Friedmann-Robertson- Walker (FRW) cosmologies filled with two ideal barotropic 
fluids with constant state parameters, where one of the energy densities is given by a sum of two powers of the scale 
factor. For the purposes of this paper, we shall define a scaling cosmology as one in which the energy densities pi and 
P2 scale exactly as powers of the scale factor in the form pi ~ a" and p2 ~ a" with n constant, leading to the ratio 
of energy densities of the form r — pi/p2 — const. On the other hand, we shall use the term quintessence not only 
for a matter content in the form of a scalar field as is generally used, but also for any kind of substance having the 
equation of state — l<a;<— 1/3. 

The outline of the present paper is as follows: In Section II we obtain the general solution for the chosen form of 
one of energy densities. In Section III we study some specific two-fluid interacting cosmologies, considered before in 
the literature, that are included in the obtained cosmological models as particular cases. In Section IV some features 
of the studied cosmological models are discussed and a three-fluid interpretation is given. Finally, Section V presents 
some concluding remarks. 



II. FIELD EQUATIONS FOR TWO FLUID COMPONENTS 



As we stated in the previous section we shall postulate the form of one of the energy densities to build a viable 
cosmological model which is able to lead to an accelerated expansion and solves or alleviates the coincidence problem. 
In order to do this the model must successfully reproduce some important cosmological scenarios considered before 
in standard cosmology. For example, FRW cosmologies fllled with a single fluid p with equation of state p = cup have 
an energy density of the form p = pQa~^^^^^\ and those with two non-interacting fluids, as we have seen above, have 
the energy densities given by ([7]). Another type of FRW cosmologies, which provides us with a way to approach the 
coincidence problem, has energy densities which scale exactly as powers of the scale factor, so they are considered to 
be given by pi = pioa" and p2 = p2Qa^ in order to have pi/p2 ^ a'^""**. Just these types of cosmological scenarios 
we want to generalize. 

The starting point is the assumption that one of energy densities, say pi, depends on the scale factor as 

pi{a) = Cia'^ + C2a^ , (9) 

where Ci, C2, a and /3 are constants. Then substituting the expression ([9]) into the conservation equation Q we 
obtain for the other energy density 

a + 3{l + uj2) p + 3(1 + ^2) 

where C is a constant of integration. It is clear that expressions (|9]) and (fTO|) may lead to an asymptotically scaling 
behavior of energy densities at early and late times, depending on which term in (a) and p, (a) is dominating at 
each epoch. 

Normalizing expressions ^ and ([T0|) to the present values pio and p20 (the present value of the scale factor is 
normalized as a = 1) we can express the constants Ci and C2 through the model parameters as 

Ci = , HP + 3(1 + ^i))pio + W + 3(1 + C.2))(P20 - C)], (11) 

^2 = ^ [(" + 3(1 + ^i))Pio + (a + 3(1 + 0.2))(P20 - C)], (12) 

3(wi - W2j(p - a) 

respectively. 

Now, by taking into account Eqs. (|9]) and (fTO|) . we can find the scale factor by solving the Friedmann equation (|3]) 
which takes the form 

k 

iH"^ + ^ = Kptotal, (13) 



where 



- n - r«-3(l+-2) , 3(^2-^1) p „a , 3(a^2 - ^l) ^ 8 

Ptotai - Pi + P2 — t-a ' -\ 7 Oifl + n , on I z^2a'^- 

a + 3(1+^2) p + 3{l+uj2) 
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Note that for C = p20 the energy density (|10l) takes the foUowing form: 

n(n\~n ^-3(1+^2) , (« + 3(1 + a;i))(/? + 3(1 + ^i))pio „ 

P2W - p2oa ' -\ TT- —T715 — -\ l« -a). (i4j 

3(a;i - W2)(p - a) 

III. INTERACTING FRW COSMOLOGIES 

As we stated above, in recent years the interacting interpretation of the conservation equation ([4]) has received 
increasing attention, which is mainly preferred by cosmofogists because it helps to study the coincidence problem by 
generating two-fluid cosmological solutions. Clearly, for the solution found in the previous section, we can introduce 
into consideration an interacting term by computing Q{t) with the help of Eq. Thus, by substituting Eq. © into 
the first expression of Eq. ([8|) (or Eq. (fTO|) into the second expression of Eq. ([8])), we obtain that the interacting term 
associated with our solution © and pHt is given by 

Q{t) = {a + 3(1 + LUi))HC\a°' + {(3 + 3(1 + uji))HC2a^. (15) 

Note that we obtain the case of two non-interacting fluids ([7]) by taking C2 — and a + 3(1 +a;i) = 0, which implies 
that the interaction term Q — 0. 

Notice that from the general solution ^ and (ITOl) we have that for C = and a — —3 (l + ^(wi + W2)) the energy 
densities take the form 

where the interaction term is given by 

Q = Co(c<^i-c^2)a-4"^^"^+"^^. 

Clearly, if wi — UJ2 we obtain the standard case of a FRW cosmology filled with a single fluid. The interaction in this 
case exists only due to the presence of different pressures pi = wip and p2 = ijJ2P- 

Let us now impose on the general expressions for energy densities ([9]) , (|10l) and the interacting term p3|) some types 
of specific interactions considered in the literature. In general the considered interacting terms are functions of the 
energy densities multiplied by a function with units of inverse of time (generally H{t) — a/ a). As we quoted in the 
introduction some often considered types are: Q — aHp^^j, Q — dtHp^^^, Q = aH{pj^^j-\-p^j^), Q ~ aH p^j^.^+PH p^^^, 
where p^^^j and are the energy density of dark matter and dark energy respectively. We will therefore consider, 
in the rest of this section, these types of interactions in order to link them associated with our solution interacting 
term in 



A. Imposing an interacting term proportional to the Hubble parameter and to one of the energy densities 

Let us first consider that the interacting term Q{t) is given by 

Q{t)^aHpu (16) 

where a is a dimensionless constant parameter. Thus, imposing on the obtained cosmological solution the condi- 
tion we obtain 

Q{t) = aHpi = aH (Cia" -I- C2a'^) ■ (17) 
Comparing Eqs. (jl5p and (jl7p we obtain the following constraints: 

a = a + 3(l+tJi),a = /3 + 3(l+wi), (18) 

implying that 

a = ^ = a- 3(H-a;i). (19) 

This implies that the energy densities now may be written in the following 5-parametric form: 

Pi(a)=Pioa"-'^'+"^\ (20) 

P2{a) = P2oa-^^'+^^^ + , ,7^° r (a-^(i+-^) - a^-'^'+^^A . (21) 

a + S[uj2 — wi) V / 
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Let us now suppose that the interacting term is given by 

Qit) = pHp2, 



(22) 



with (3 a dimensionless constant parameter. Thus, by taking into account Eqs. (fTO)) . (|T5|) and (|22|) . we obtain the 
foUowing condition: 



(a + 3(1 + Wi))Cia" + (/3 + 3(1 + uj{))C2af^ = /3 I Ca 



a + 3(l +a;2) 



/3 + 3(i + c^2) y ^ ^ 



It is clear that we must put C = in order to have a self-consistent interacting solution and the foUowing constraints 
must be imposed: 



implying that 



Thus from Eqs. 



^ = -a - 3(1 + W2), /3 = -/3 - 3(1 + W2), 

a = /3 = -/3- 3(1+^2). 
and (|10p we conclude that these interacting scenarios may be written as 



Pi(a)=pioa-^-3(i+-), 
P2(«)= ^^"^"J^^"^ Pio«^W^.)^ 

where pio = Ci + C2. It is clear that in this case we have a scaling behavior of energy densities where 

_ P2 3(lUi ~ UJ2) ~ P 



Pi 



(24) 
(25) 

(26) 
(27) 

(28) 



This ratio is positive for < /3 < 3(a;i — ^2) when wi > 102, and 3(aJi — ^2) < /3 < when uji < lj2. For /3 = 3(a;i —UJ2) 
we obtain the standard FRW solution for a single ideal, and for ui — lj2 we have that pi = —p2^ thus obtaining a 
vacuum FRW cosmology. 

Nevertheless, notice that the 4-parametric solution (|26p and (P7)) is a particular case of Eqs. ([201) and ([?!]) . Effec- 



tively, Eqs. (PSI) and (P7)) are obtained if we put p20 



and + 3(^1 - ^2) into Eqs. (HOI) and (gl 



Therefore, since the coupling (1^^ is a particular case of Eq. ([TS]). the general solution for a coupling proportional to 
the Hubble parameter and to one of the energy densities is given by Eqs. (|20|) and ([2T|) . 

Now we shall consider some particular cases of interacting cosmological scenarios described by Eqs. ([20)1 
and (12111 in order to confront them with those discussed in the literature. 



1. Interacting term proportional to dark matter energy density 

Let us first suppose that pi describes the energy density of the dark matter component p^^^ . This implies that 
we must put wi — 0. In this case the interacting term (jl6p may be written as Q{t) = aHp^^^j and the energy 
densities (PU)) and (PT|) as follows: 

PnM ■= Pi (a) = PnMo a""^ (29) 

Pde ■- P^W - Pdeo'^ + 3^2 V /' ^ 

where Pj-,^j„ and Pj^^,,, are positive constants. For —1 < uj2 < ~l/3 the dark energy p2 may be interpreted as 
quintessence, and for ^2 < — 1 as phantom matter. Clearly, due to the interaction (|16p . Eq. (|29p is a deviation from 
the standard behavior of the matter component Pj^^^ ^ which implies that matter is conserved separately from 
P2 ■ Some aspects of this kind of interacting scenarios were considered and discussed by authors of Ref. 0] ■ Let us 
now enumerate other relevant features of these cosmologies: 

By supposing that the dark matter energy density decreases with the expansion we see from Eq. (|29p that we must 
require that a — 3 < 0. On the other hand, from now on we will be considering only positive energy densities. Thus, 
in order for the energy densities (|29p and (j30l) to correspond to positive densities during all evolution, we must require 
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that - ,^ < and Pn„n + "''F^° > 0. The condition -:r-^ — < imphes that a > and a + 3oj2 < or that 5 < 
and a + 3a;2 > 0. Any other combination for a and UJ2 will necessarily imply that the energy density of dark energy 
becomes negative for some value of the scale factor a during the cosmic evolution. 

Note that if a < and a + 3a;2 > we have that aj2 > and then in this case the consideration of dark energy is 
excluded. There exists the possibility of considering uj2 < ~l/3 for constraints a > and a + 3aj2 < which imply 
that a;2 < 0. In this case the coupling parameter a is constrained to be < a < 3, implying that in these scenarios 
the energy is transferred from dark energy to dark matter. At early times we have that the dark energy behaves as 

P2 — — ^'^'i then r = PdmI Pde — — (a + 3a;2)/Q!, and for late times r — > implying that the dark energy 

dominates over dark matter. For this model one finds that the equilibrium between dark matter and dark energy, i.e. 
r(aeg) = 1, corresponds to 

_ f ^ g + 3cj2 PoEo \ 

""^ V2a + 3w2 2a + 3w2PDMo/ 

while the accelerated expansion begins at 

_/ a(i + 3c^2) ^ (i + 3^2)(a + 3cj2) P,,^o \^^ 

\'iuj2(a-\) 3w2(a-l) Pdmo) 

Clearly these two expressions do not coincide. It can be shown that for W2 < —2/3 we have a^c < a-eq implying that 
we can have scenarios where the universe already has entered into an accelerated expansion while the dark matter 
component is still dominating. For —2/3 < 0^2 < —1/3 we have that aac > deq and then the acceleration begins when 
the dark energy component already dominates in the Universe. It is interesting to remark that this kind of behavior 
also is observed in power law cosmologies studied in Ref. fs"] 

Let us note that, in this case, in order to have a transfer of energy from the dark matter component to dark energy 
we need to require d < 0, then having a negative dark energy at early times. 



2. Interacting term proportional to dark energy density 

It is interesting to note that the considered energy densities (|20|) and pip are not symmetric for the state parameter 
interchange {uji = 0,a;2 = i^) ^ ('^i = uj,uj2 =0), where w is a free parameter, as well as for the general solution Q 
and pop . So for scenarios containing interacting dark energy and dark matter with cj2 = the energy densities are 
given by 

:=pi(a) = p„,„a"-^(i+"^), (31) 
Pn. P2ia) = P^.oa-^ + 1^ - „^-3(i+..)) . (32) 

The interpretation of this qualitatively different interacting cosmology is direct: p2 is now the energy density of dark 
matter while pi is the energy density of dark energy for uji < —1/3. Thus from Eq. (|16l) we conclude that now the 
interacting term is proportional to dark energy density, i.e. Q{t) = aHp^j^. This solution is a particular solution of 
the interacting scenarios discussed in Ref. Q (see Sec. IIA). 

Now we shall enumerate some relevant features of these cosmologies. Note that in order to have a positive energy 
density for dark matter during all evolution we must require that (5/ (a — 3a;i) < and p^j^j„ + {a/ {a — 3uji))p^j^„ > 0. 
The condition a/{a — 3uJi) < implies that a > and a — 3lui < or that a < and a — 3uji > 0. On the other 
hand, by supposing that the energy density of dark matter (1321) decreases with the scale factor we must require that 
a-3(l+wi) <0. 

For a > and a — 3aji < we conclude that wi > and then the consideration of dark energy is excluded 
in this case. On the other hand, for a < and a — 3wi > we conclude that wi < 0, implying finally that 
the constraint on the coupling parameter a must be 3a;i < d < 3(1 + uji). In this case at early times the dark 
matter component dominates over dark energy and for late times we have a scaling behavior for energy densities since 
r — Pom I Pde — —0!/{d — 3loi). In such scenarios the transfer of energy goes from the dark energy to the dark matter 
component. 

In general for these models we have that the equilibrium between dark matter and dark energy corresponds to 

a ^ ( " ^ 3cji - a Po^„, \ 
\2d- 3a;i 3a;i - 2d p^Eo ) 
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while the accelerated expansion begins at 

1 

^ _ f " I (a -3^1) PoMo Y'""^' 

\3a;i(l + Scji - a) 3wi(l + 3wi — a) p^^.^ / 

As before we can have scenarios where the universe already has entered into an accelerated expansion while the dark 
matter component is still dominating. 

Lastly, notice that in this case in order to have a transfer of energy from dark matter to dark energy we must require 
a negative dark energy during all cosmic evolution, or a > with a negative dark matter at early or late times. 



B. Imposing an interacting term proportional to the Hubble parameter and to a linear combination of the 

energy densities 

Let us now consider the interacting term given by 

Q{t)^&Hpi+pHp2. (33) 

Thus from Eqs. ©, PH)) . ([T5|) and ([55)1 we firstly conclude that we must put C = 0. Thus the constraints on the 
model parameters are 

^= . on , ^+a + 3 1+ wi ,a= o.n,^, ^ + /? + 3 1 + wi . 34 



Note that these conditions include the cases ([16]) and ([22|) considered above since for a = we obtain the condi- 
tions ([Ml) , and for /3 = we obtain the conditions ([T5[) . However note that in the procedure below we assume that 
/? ^ in order to self-consistently solve the now treated problem. 

In the following we shall express the powers of the scale factor a and (3 as functions of the model parameters wi , 
a;2, a and (3. Thus from Eqs. ([M|) we conclude that the powers of the scale factors are constrained to be given by the 
expressions 



" =-P+ =-P+ ^e,, (35) 



where = ±1, = ±1 and 



g = (a - /3) + 9 (cji - cja)' + 6 (a + /?) (wj -uji),p=-- (5 - /s) - - (wi + + 2) . (36) 
Thus the interacting cosmological scenarios take the following forms: 

pi{a)^CiaP+-2^^V^ + C2aP+'^'fi^ (37) 

P2{a) = _P±KV^+^il±^ Ca^^h^^^^- P+h.y-^ + ^i^^-^) c^a^H^.^., (38) 
P+5e„y9 + 3(l+W2) p+ie,Vg + 3(l+W2) 

Clearly we can have scenarios where a 7^ /3 by taking = ^^0' ^^"^ scenarios where a = /? by taking = or 

by requiring (7 = 0. This latter condition will imply that a — 3(a;i — UJ2) + /3 ± 2-^3 P{uJi ~ ^2) and then the energy 
densities take the following particular form: 



pi(a) ^ i^a-3(i+-2)±V3^("i— 2)^ 



£> ^3(^1 - ^2) ± iP{uJi ~ LJJ2) 



P2{a) = T— ^ ^ X a-3(i+-2)±V3^(-i— 2)^ (39) 

y3/3(wi - W2) 

where is a new constant. 

This kind of solutions was considered by authors of Ref. 10]. Our solution takes the form of the solution discussed 
in Ref. [lo| by means of Wi = — 1, ^2 = 7b — 1, 5 = — /?b and f3 = a^, where all parameters with subscript 
B denote the parameters used by Barrow and Clifton in cited references. However, it must be noticed that in this 
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case we can not have a scaling behavior for the energy densities pi and p2 of the form p2 = X{C3a^ + C4a^) = Xpi, 
with A, C3, C4, u and v constants. This can be seen by direct integration of Eqs. ([5]) with the interacting term p3p 
by imposing the condition p2 — Xpi with A constant. In this case the general self-consistent solution is given by 
Pi — piqu ^ 1+^ '. This scaling solution also can be directly obtained from the general solution (|9]) and (|10]) by 
imposing the same condition p2 = Xpi- 

Another scenario considered in interacting cosmologies is that defined by an interacting term of the form Q{t) = 
aH{px + P2). These cosmologies are a particular case of Eqs. (|33)) and (p4)) and may be obtained from Eqs. (|35|) - (|38)) 
by putting a = /3. In general this type of interacting terms has been discussed in the framework of scalar fields [11]. 

Lastly notice that the solution ((37)) and (|38| does not include as a particular case the whole class of interacting 
cosmologies given by and described by Eqs. (^0]) and (^1]) . This is due to the fact that for /3 = we must impose 
the extra conditions (3 + 3{1 + UJ2) — 1 and a + 3(1 + UJ2) = 1 as we can see from Eqs. (p4|) . thus obtaining only 
particular solutions of Eqs. ([20]) and ([21]) . For the case 5 = Eqs. ([34]) do not impose any extra condition on the 
model parameters. Thus the whole class of interacting cosmologies (1^^ and ([^T)) is included in Eqs. (|57)) and psp as 
a particular case. 



IV. SOME FEATURES OF THE GENERAL SOLUTION 



Now we shall study the cosmological scenarios derived in Sec. II by the assumptions that one component of the 
universe obeys an equation of the form ([9]) and that pressures obey barotropic equations of state with constant state 
parameters. As we have seen cosmologies described by Eqs. ^ and (fTU]) . with arbitrary a and /3, generalize the 
previous interacting cosmological scenarios considered in Sec. III. So clearly the found analytical cosmological models 
provide us with extensions of the analysis of cosmologies containing dark matter and dark energy. 



A. Cosmologies with cji = 



Let us first consider generalizations of the model described by Eqs. (P^)) and ([50)) . This implies that we must put 
cji = into Eqs. (O and Note that for C2 = 0, Ci = p^.^, C = + |^f|^ and a = d - 3 (or Ci = 0, 

C2 = Pdmo : = PoEo + ^^^^ /? = a — 3) we obtain the same cosmologies described by Eqs. ([29]) and ([30]) . Thus 

the extension of the analysis is obtained for Ci ^ and C2 ^ simultaneously. In this case Eq. ^ describes the 
behavior of the dark matter and Eq. ([10|) describes the behavior of the dark energy. By taking into account Eqs. ([TT|) 
and (|12p (with pig = Pdmo ^"^"^ P20 = Pdeo ) dark energy component now is given by 



{a + i){fi + i)p _ ^ p - C ^ ^ ^ „ + 3)(^ + 3(1 + . (40) 

3a;2(a-p) 3w2(a-p) 

Clearly, in general, there are values of the model parameters which lead to negative values of the dark matter 
and/or the dark energy component. These energies can be either positive during all cosmic evolution or negative at 
the beginning or at the end of the expansion of the Universe. It must be noticed that negative values of dark energy 
could in principle be allowed if, for example, this type of energy is a manifestation of modified gravity [12] . 

In the following for the sake of simplicity we shall consider the case C = Pdeo- Clearly this particular form of 
the general solution is still a generalization of Eqs. ([29)) and ([30]) . which can be obtained by putting a = a — 3 and 
P = —3(1 -1-^2 ). By requiring the positivity of the dark matter the constraints Ci > and C2 > must be fulfilled. 
On the other hand, in order for dark matter to decrease during all cosmic evolution we must require that a < and 
/3 < and, since 012 is the state parameter of the dark energy, we shall constrain it to uj2 < —1/3. Without any loss 
of generality we can consider a > (3. Thus we have that p^^^^ > during all cosmic evolution if i) —2 < /? < a < 0, 
<UJ2< and -2 < a < 0, -3 < /3 < -2 (in this case < 0J2 < -1/3) or ii) -3 < a < and l3 < ~3 

for any lu2 < —1/3. 

Now, by taking into account Eq. ()15|) . we conclude that for the case i) we have energy transfer from dark energy 
to dark matter (Q > 0), while for the case ii) we have that at early times the energy is being transferred from dark 
matter to dark energy (Q < 0) and at late times we have a transfer of energy from dark energy to dark matter 

(Q>o). ^ ^ 

Lastly notice that, in general, independent of which term among a", and dominating at early 

and at late times, we may have a scaling behavior for energy densities at early and late epochs. However, we are 
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interested in such a behavior of the considered set of solutions at late times. Without any lost of generality we can 
suppose that /3 < a < 0. Thus, in this case the dark matter component will behave as Pj-,j^j ^ Cia". For the dark 
energy component we can have p^^ ^ Ca^'^(^+"^' for a < —3(1 + W2), or ^ Ca" for —3(1 + LJ2) < a, where 
C'^ 34i^[(/5 + 3KMo + (/? + 3(1+^2)) (P^.o -C^)]- Thus, for a < + we have r^^^ ^a"+3(i+-.). 
Clearly in this case r vanishes at late times, implying that the dark energy dominates over the dark matter. For the 
second case we have that 

^ _ Pdm ^ q; + 3(1 +1:^2) _ _-|^ _ 3i:j2 
y0j3^ a + 3 a + 3 ' 

thus alleviating the coincidence problem. In this case, in order to have r > we must require —3 < a < 0. 



B. Cosmologies with tJ2 = 



Now we shall consider the generalization of scenarios ([31]) and ([32|) . This implies that we must put W2 = and 
then Eq. (jlOp describes the behavior of the dark matter component, while the energy density of the dark energy p^^ 
is described by Eq. ([9]), with uji < —1/3. Thus the dark matter component has the following form: 

_ ^ -3 ^ (t^ + 3(l + a;i))(/3 + 3(l + ^i)K^„ . „ 



3ciJi(/3 — a) 

Notice that in general, as in the previous case, there are values of the model parameters which lead to negative values 
of the dark matter and/or the dark energy component (cither at early or at late times). 

In the following for the sake of simplicity we shall consider the case C = Pdmo- clear that this particular form 
of Eqs. ^ and (HT]) contains as a particular case the solution (|5T|) and ([5^ . This can be seen directly by taking 
a = a — 3(1 + wi) and (3 = —3 (this condition leads to C2 = 0, Ci = Poeo)- order for dark matter to drop during 
the expansion we shall require that a < and (3 < 0. It is clear that for wi < —1/3 there is a set of values of the 
model parameters which leads to a positive energy density of dark matter during all cosmic evolution. It can be 
shown that for a < —3 and /? < —3 the energy density of the dark matter takes negative values at early times of the 
evolution. We may find a set of values for the parameters which leads to a positive dark energy during all evolution 
if —3 < a < and —3 < /? < 0. For example in Fig. [l]we show a case where the energy density of dark matter takes 
only positive values during all evolution. It is also shown the behavior of the interacting term Q, which is positive at 
the beginning of the expansion, becoming negative at some value of the scale factor. Note that in Fig. [T] is considered 
a flat FRW scenario (fc = and then rij„j„, = 1) with fi^jj^o = 0.3. This value of the dimensionless energy density of 
the dark matter is in good agreement with a wide range of observations: high-redshift Type la supernovae, evolution 
of galactic clusters, high baryon content of clusters, lensing arcs in clusters, and dynamical estimates from infrared 
galaxy surveys |13l - [la | . It is interesting to note that the interacting term ([TS]) now is given by 

^^^^ ^ (a + 3(l + ^,))(^ + 3(l + u.,))p.o [_(^ ^ 3)^.-1 ^ ^ 3)„,-i] (42) 

Thus we have that it does not change sign only if a < —3, /3 > —3 or a > —3, /3 < —3. However in this case there 
are configurations leading to a negative energy density of the dark matter component. 

It is interesting to remark that, as in Sec. IV-A, for any solution with W2 = we may find a set of cosmologies 
presenting a scaling behavior for energy densities at early and late epochs. Effectively, without any loss of generality 
let us suppose that /3 < a < 0, and then at late times the dark energy component will behave as p^j^ ^ Cia" while 

the dark matter component as p^j,^ ^ Ca~^ for a < —3, or pj^j,^ ^ Ca" for a > —3, where C — ^5^7p^^^[(/3 + 3(1 + 

^i))PnEo + {PnMo ~ + 3)]- Thus for a < — 3 we have for the ratio of energy densities r = ~ ■^a~''"~^^\ 

Pde ^1 

Clearly in this case r diverges at late times, implying that the dark matter component dominates over the dark energy, 
so these models must be ruled out of consideration. On the other hand, for a > —3 we have cosmological models with 
a scaling behavior since 

Pom a + 3(1+^1) 3wi 

J" = = — 1 — 

PoE a + 3 a + 3 ' 



thus alleviating the coincidence problem. In this case, in order to have r > we must require a + 3(1 + ui) < 0. 
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FIG. 1: We show the evolution of the dimensionless energy densities Hi for dark matter (solid line), dark energy (dotted line) 
and the corresponding behavior of the coupling term Q (dashed line) vs. redshift z for the following values of the model 
parameters: k = 0, a = —2, /3 = —1, Hj^^^q = 0.3, H^eo ~ O-^i ~ —0.9. In this case the interacting term Q is negative for 
z < 1, implying that for z > 1 we have a transfer of energy from the dark matter to the dark energy, and for z < 1 the energy 
is being transferred from the dark energy to the dark matter component. 



C. The three-fluid interpretation 

We next consider models in which there are two dark matter sectors: one describing the standard visible matter 
which is not interacting with the dark energy component, and another which is interacting with the dark energy. Both 
of these dark matter components are treated as pressureless fluids. The advantage of such a description is that we 
can describe the behavior of visible matter with an energy density which decreases with the scale factor as p^n a , 
while the energy density of the another dark matter component behaves as p^^, ^ a", with n < and n ^ —3, in 
order that this second dark matter component also decreases with the scale factor. 

Let us consider the general solution ^ and pUj) with uj2 = 0. In this case we shall not consider the energy density 
P2 in the normalized form (|4T|) and we shall write it in the following form: 



P2^PuM+Pm, (43) 



where 



a + 3(l+c^ a ^ + 3(1+^ p 
Pdm = -T. -^—7. , 

Prn = PmOa"^, (44) 

and we have put C — pmo > 0. Note that, as in Sec. IV-B, the energy density of the dark energy is given by := pi 
of Eq. ©. 

Thus, by putting expressions (|44| and (|9]) into the conservation equation (|4]), it may be rewritten as 

Pdm + Pde + ^H{p^^, + (1 + ^oe)Pde) ^ 0, 

Pm + iHp,n = 0, (45) 

where we have put ui := w^^, . Clearly the interpretation of these equations is different here: we have the conservation 
of the total visible matter pm, while the dark matter is conserved together with the dark energy component. Therefore 
we can associate an interaction between p^j^,^ and p^j^ . In this case the coupling between dark matter and dark energy 
is still given by Eq. (fT5|) . 
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FIG. 2: We show the behavior of the dimensionless energy densities Q,i for standard visible matter (dotted hne), dark matter 
(dashed hne) and dark energy (solid line) for the following values of the model parameters: fc = 0, QmO = 0.04, J^oa/o ~ 0.26, 
oil = -2/3. This implies that /3 < -1.542 and -1.542 < a < -1. We have plotted curves for a = -3/2 and /3 = -5/2. For 
late times the dark energy component dominates, while at early times dust is dominating. 



Normalizing the solution in order to have today p^^.j{a = 1) = Pj^^jg and Pj^^i^ ~ 1) ~ Pdeo the constants Ci and 
C2 niay be written as 



C2 = 



{a + 3) [(/3 + 3)p„,,„ + (/? + 3(1 + c^i))p„, 

3wi(a-/3) 
{(3 + 3) [{a + 3)p^,,„ + (q + 3(1 + 

3cJi(a — 13) 



(46) 



It is clear that, if we want to have positive energy densities during all cosmic evolution, we must require that all terms 
of Eqs. ^ and ^ must be positive. Thus, by requiring Ci > 0, C2 > 0, < 0, < Q and without 

any lost of generality /3 < a we have that the following constraints must be satisfied for negative values of the state 
parameter ui. 



/3 + 3(l+LJi) ^ p^ 



< 



Pdmo ^ a + 3{1 + UJi) 



< 



a + 3 



< 0. 



This implies that the parameters a and /3 are constrained to be 

3(l + ro+wi) 



1 + ro 

-3 < /? < 



< a < -3(1 + wi), 
3(1 + ro +a;i) 



1 + ro 



(47) 



where rg = ""^ . It is clear that these constraints imply that energy is being transferred from dark energy to dark 

^DEO 

matter. Notice that for Ci = or C2 = we have that the energy densities of dark matter and dark energy are 
functions of the same power of the scale factor, i.e. PoeI Pdm — const. 

In Figs. [5] and 13] are plotted the evolution of three dimensionless energy densities i7„i, Vt^j^, and fi^E {^m + ^da/ + 
= 1) in the case of flat FRW cosmologies {k = 0) and always positive dark energy and dark matter energy 
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FIG. 3: We show the behavior of the dimensionless energy densities Q,i for standard visible matter (dotted hne), dark matter 
(dashed hne) and dark energy (solid line) for the following values of the model parameters: fc = 0, QmO = 0.04, J^oa/o ~ 0.26, 
coi = —1.2. This implies that /? < —0.375 and —0.375 < a < and we have plotted curves for a = —0.3 and /3 — —1/2. For 
late times the dark energy component dominates, while at early times dust is dominating. Note that in this case the dust-dark 
energy equality and dust-dark matter equality occur for z 2 and z ^ 1 respectively. 



densities. In order to have an estimation for the behavior of dimensionless energy densities, the cosmological model 
parameters are fixed to take the values fJ^^ = 0.04 and = 0.26, which are in good agreement with a wide range 

of observations [l3l - [l5| . In both cases the energy is transferred from the dark energy to the dark matter. 

If we do not exclude the possibility of negative energies we can also consider cases where the dark energy becomes 
negative at some stage during the cosmic evolution. For example in Fig. |4] we show the case where negative values 
for the dark energy density are allowed at early times. However it must be noticed that the total energy density 
Ptot ~ Pmo + Pdmo + Pdeo positive during all cosmic evolution. An interesting feature of such a scenario with 
negative energy is that we have a period where the energy is transferred from dark energy to dark matter and another 
where this occurs from the dark matter component to dark energy. 



V. CONCLUSIONS 



In this paper we have studied Friedmann-Robertson- Walker cosmological models with matter content composed of 
two barotropic perfect fluids where one of the energy densities is given by a sum of two powers of the scale factor of 
the form of Eq. By associating with these cosmologies an interacting term Q, it can be shown that interacting 
scenarios with couplings given by Q = aHpi, Q = aHp2, Q = aH{pi + P2) and Q = aHpi + PHp2 (with constants 
a and /?) correspond to particular cases of our cosmological model. The studied cosmological models contain a class 
of solutions having a scaling behavior at early and at late times, and then the coincidence problem is substantially 
alleviated. 

It is interesting to note that in the framework of the considered scenarios it is possible to introduce a three fluid 
interpretation, where one fluid describes the standard visible matter and is conserved separately from the dark matter 
and dark energy components, which are conserved together. If we want to introduce the interacting picture, thus 
the visible matter is not interacting neither with the dark matter nor dark energy components, while these latter 
two components are interacting with each other. In this case, if energy densities of dark matter and dark energy are 
positive during all cosmic evolution, then the energy always is transferred from the dark energy component to dark 
matter in agreement with the conclusions of Ref . [16| . 

It is remarkable that the inequalities (j47p impose a lower limit on the values of the state parameter wi . Effectively, 
if we require that the energy density of the dark matter be decreasing with the expansion, i.e. a < and /? < 0, 
then the state parameter is constrained to be in the range wi > — (1 + j'q). Thus, in agreement with the observations, 
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FIG. 4: We show the behavior of the dimensionless energy densities Q,i for standard visible matter (dotted hne), dark matter 
(dashed line) and dark energy (solid line) for k = 0, QmO = 0.04, ^j^^^g = 0.26, toi = —1.2, a — —0.275 and P = —3.1. In this 
case the dark energy density becomes negative for z >9 and the dust-dark energy equality occurs for z ~ 1.5, while for z ~ 2 
we have the dark energy-dark matter equality. Notice that for z > 5 the energy is transferred from the dark matter to dark 
energy, and for z; < 5 is from dark energy to dark matter. 

we have for the ratio tq — ^smsl — ^dmo _ 0.26/0.7 and then uji > —1.371. Note that this constraint on the state 
parameter uii is allowed in this model by the requirement of positivity of energy densities of dark matter and dark 
energy during all cosmic evolution. However, one can impose more stringent constraints on the model parameters by 
using for example the latest supernova data. In order to do this we shall use the recent compilation of SNe la data 
called the Union data set [13, which contains a set of 57 nearby (0.015 < z < 0.15) Type la supernovae, and 250 
high-redshift supernovae. The statistic is quite helpful in constraining the parameter values of a given model [l^. 
In our case this method will allow us to fit the set of the model parameters a, (3 and wi to the set of cosmological 
parameters of the Union compilation data, by finding the best fit values of the model parameters by minimizing 
this x^- In this case the minimum of should be roughly equal to the number of data, or the so-called "reduced 
chisquare" xl = Xmin/'^^'^ should be roughly equal to 1. 

In Fig. [5] we show the probability contours from the above 307 data points only at 68.3%, 95.4% and 99.7% 
confidence levels (from inside to outside) in the a — uii plane for the value f3 = —5/2 of Fig. [5] In this case the 
best-fitting parameters are a — —0.1 and uji — —0.97 with xt — 1.021, implying the constraints —3.6 < a < 3.4 and 
— 1.49 ^ < —0.45 (99.7% C.L.). Note that in this case, by taking into account the above constraint on the state 
parameter ioi, we can rewrite the latter constraint as —1.371 < wi < —0.45. 

In Fig.|n]we show the probability contours at 68.3%, 95.4% and 99.7% confidence levels (from inside to outside) in 
the a — (3 plane for the value wi — —1.2 of Fig. [3] and Fig. 21 In this case the best-fitting parameters are a ~ 0.59 
and /3 = -2.99 with xl = 1-024, implying the constraints -0.34 <a< 1.52 and -4.65 < ^ < -1.33 (99.7% C.L.). 
Thus, in this case, in order to have a decreasing with expansion dark matter energy density, we must require that 
-3 < /3 < -1-33 and -0.34 < a < 0. 

Lastly, in Fig. [7]we show the evolution of the Hubble parameter H{z)/Hq vs. the redshift z for coi — —1, while in 
Figs, m |9] and [10] we show the evolution of the deceleration parameter q vs. redshift z for some values of the state 
parameter wi (—1.371 < wi < —1/3). Note that in Fig. [TUl for comparison, the prediction of the ACDM model is 
also shown. 
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FIG. 5: We show the probabiUty contours at 68.3%, 95.4% and 99.7% confidence levels in the a — wi plane for the value 
P = —5/2 of Fig. [2] In this case the best-fitting parameters are a = —0.1 and cji = —0.97 with xi = 1.021. 
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FIG. 6: We show the probability contours at 68.3%, 95.4% and 99.7% confidence levels (from inside to outside) in the a — P 
plane for the value uji = —1.2 of Fig. [2l In this case the best-fitting parameters are a = 0.59 and /3 = —2.99 with xi ~ 1-024. 
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FIG. 8: We show the evolution of the deceleration parameter q with the rcdsliift z for k = Q, f^mo = 0.04, ^^mo ~ 0.26, 
a = —3/2 and p = —5/2, and oji = —0.55 (dashdotted line), ui\ = —0.58 (dotted line), wi = —0.6 (solid line) and oji = —2/3 
(dashed line). In this case the value of the redshift z, for transition from decelerated expansion to accelerated expansion, 
becomes smaller with the increase of the state parameter wi. 
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FIG. 9: We show the evolution of the deceleration parameter q with the redshift 2 for fc = 0, Slmo — 0.04, 0,^^^^ = 0.26, 
/3 = —5/2, uj\ — —2/3, and a = —3/2 (dashed line), a — —1.2 (solid line), a — —1 (dotted line) and a — —0.9 (dashdotted 
line). In this case the value of the redshift z, for transition from decelerated expansion to accelerated expansion, becomes 
smaller with the increase of the parameter a. 
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FIG. 10: We show the evolution of the deceleration parameter q with the redshift z for k = Q, Q.mo = 0.04, ^0^0 = 0.26, 
a = —3.1, /3 = —0.05, and wi = —1.05 (dashed line), oji = —1 (dotted line) and wi = —0.95 (dashdotted line). For comparison, 
the prediction of the ACDM model (solid line) is also shown. 



